Spin-dominated waveforms for unequal mass compact binaries 
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We derive spin-dominated waveforms (SDW) for binary systems composed of spinning black holes 
with unequal masses (less than 1 : 30). Such systems could be formed by an astrophysical black hole 
with a smaller black hole or a neutron star companion; and typically arise for supermassive black 
hole encounters. SDW characterize the last stages of the inspiral, when the larger spin dominates 
over the orbital angular momentum (while the spin of the smaller companion can be neglected). 
They emerge as a double expansion in the post-Newtonian parameter e and the ratio £ of the orbital 
angular momentum and dominant spin. The SDW amplitudes are presented to (e 3//2 , £) orders, while 
the phase of the gravitational waves to (e 2 ,£) orders (omitting the highest order mixed terms). To 
this accuracy the amplitude includes the (leading order) spin-orbit contributions, while the phase 
the (leading order) spin-orbit, self-spin and mass quadrupole-monopole contributions. While the 
SDW hold for any mass ratio smaller than 1 : 30, lower bounds for the mass ratios are derived 
from the best sensitivity frequency range expected for Advanced LIGO (giving 1 : 140), the Einstein 
Telescope (7 x 10~ 4 ), the LAGRANGE (7 x 10~ r ) and LISA missions (7 x 10~ 9 ), respectively. 

PACS numbers: 04.25.Nx, 04.30.Db, 04.80.Nn, 95.85.Sz 



I. INTRODUCTION 

The orbital evolution of the binary system forming 
compact objects is accompanied by emission of gravita- 
tional waves. The slight dissipativc effect of the gravita- 
tional radiation allows for a long inspiral phase, followed 
by a short merger /plunge [l[ and the relaxation of the 
newly formed object in a process called ringdown 
The inspiral dynamics can be accurately described by a 
post-Newtonian (PN) expansion in terms of a small (but 
ever increasing) parameter e = Gm/c 2 r w v 2 /c 2 (with 
m = mi + 77i2 the total mass, r the orbital separation 
and v the orbital velocity of the reduced mass particle 
/i = m\railra). 

The inspiral dynamics including the contributions of 
the spins and quadrupole moments was extensively dis- 
cussed in Refs. 0-[H|- The gravitational waveforms h+ 
and h x characterizing the inspiral on circular orbits were 
given to 1 PN order in Rcf. 6] , including the leading or- 



der spin-orbit contributions. In Ref. [14j the spinning 
waveforms were given to 1.5 PN accuracy, and specified 
for the equal mass case by an expansion in the small an- 
gle i (which will be denoted a in our formalism, and will 
not be small) span by the orbital angular momentum L jy 
and total angular momentum J. Spinning gravitational 
waveforms on generic orbits were also given in Ref. [22| 
to 1.5 PN order accuracy, and corrected in Ref. [23|. The 
multipole moments including spin effects to 2.5 PN order 
for gravitational wave amplitudes were calculated in Ref. 

The mass range of neutron stars is relatively narrow 
due to the upper bound on their masses represented by 
the Oppenheimer-Volkoff limit [25[, [26} . Consequently 
the equal or nearly equal mass compact binary model 
stands as a good approximation during their inspiral. 

The situation is however radically different for black 



hole binaries, and these are the compact objects, which 
have significant spin. There is no reason to believe that 
for astrophysical black hole binaries (with each black hole 
having a mass extending from a few solar masses, M 
to a few ten times of M Q ) the comparable mass case is 
more likely than any other mass ratio v — mil 'm\. More- 
over, for supermassive black hole binaries (with masses 
3 x 10 6 — 3 x 10 9 M©) it has been shown that a mass ra- 
tio between 0.3 and 0.03 is typical @, [l(|, [27j equal 
mass and extreme mass ratio encounters being disfa- 
vored. Low mass ratio and intermediate mass black hole 
binaries emitting gravitational waves are also believed 
to be sources for the LIGO detectors. However no sig- 
nals were detected yet [H}, [H|, [H .With the advent 
of advanced LIGO 1311 . the planned space-born missions 
LISA-eLISA/NGO lil-dl and LAGRANGE [H, also 
the third generation detectors, like the Einstein Tele- 
scope [35}, the knowledge of waveforms from binaries 
with smaller mass ratio ranges becomes imperative. The 
small mass ratio parameter in these cases stands as a 
second small parameter, which modifies (and simplifies) 
some of the immediate results standing for the equal mass 
case. Indeed, it has been shown in Ref. [9|, that the ratio 
of the spin magnitudes is 

S i = ^ , (1) 
£1 Xi 

where Xi G [0, 1] are the dimensionless spin parameters. 
For rapidly spinning compact binaries and for small mass 
ratio the role of the second spin becomes negligible. 

Moreover, the ratio of the dominant spin and magni- 
tude of Newtonian orbital angular momentum is @ 

Si A/2 -1 (2) 



This relation shows that while for equal masses the in- 
spiral is dominated by Ljv , for small mass ratios this can 



2 



change. Indeed, Si becomes dominant in the last stages 
of the inspiral for the unequal mass case v < 0.1. (Obvi- 
ously, in the test particle limit Si dominates throughout 
the inspiral, the orbital angular momentum of the test 
particle being negligible.) 

In this paper we will focus on the situation when v 
is small and we will derive the corresponding spinning 
waveforms to accuracy of order v. This immediately im- 
plies to disregard S2 as of order v 2 compared to Si. 1 As 
the gravitational wave amplitude increases drastically to- 
ward the end of the inspiral, we will concentrate on this 
regime. As has been shown in Q, [l(|, in this mass ratio 
range the end of the inspiral is characterized by Si domi- 
nating over Ln- Hence we call the corresponding gravita- 
tional waveforms as spin- dominated waveforms (SDW). 
We express this dominance by introducing a new small 
parameter 

i = e-^v, (3) 

assuming £ < 0.1. As the PN parameter increases during 
the inspiral, this condition selects the last part of the 
inspiral starting from a given radius n . Hence the SDW 
are valid from r\ to the end of the validity of the PN 
expansion. 2 

We derive the SDW inspiral waveforms for circular or- 
bits as a double expansion in the parameters e and £, 
with a linear accuracy in £ and 3/2 orders in e (dropping 
however the £ 3 / 2 £ terms). The introduction of the pa- 
rameter £ leads to the natural neglection of the second 
spin, since Si — G/cmfxi such that from Eq. (fTJ) we 
find S2 = x 2 G/cm\v 2 = \ 2 G ' j 'cm\e^ 2 . Since the first 
terms containing spins in the amplitude are the 1 PN 
order terms the leading order terms containing S2 
would be shifted to 2PN orders in the amplitude, to be 
neglected in our approach (unless £ 3> 1, which is not 
what we assume in this paper). 

The structure of the paper is as follows. In Sec. ITT1 we 
analyze the limits of validity of the proposed SDW. In 
Sec. IIIII and Appendix |A"1 we derive and present the SDW 
(gravitational wave amplitudes, with the spin-orbit cor- 
rections included) resulting from the double expansion. 
Here we employ that in the discussed parameter region 
the angle j3 1 span by the dominant spin and total angular 
momentum is small, of order £. Therefore wherever pos- 
sible, we express the angle a span by the total angular 
momentum and orbital angular momentum in terms of j3 1 
and perform the respective expansions to first order in £. 
Section IIVI contains the expressions of the gravitational 



1 Gravitational waveforms with the inclusion of only one spin were 
previously investigated in Refs. [3(J, [37| |. to leading order in 
amplitude and 3.5 PN orders in the phase (the physical template 
family), however omitting the qua drup ole-monopole and self-spin 
contributions, given later in Ref. |38|l . 

2 The mass ratio v will be replaced by e 1 ' 2 ^ in all expressions, 
with the exception of leading order contributions, where we keep 
v. 



wave phase, expanded in terms of both e and £, with 
a linear accuracy in £ and second order in e (dropping 
however the £ 2 £ terms). The phase then contains spin ef- 
fects (the leading order spin-orbit and self-spin) and mass 
quadrupole-monopole effects, however spin-spin terms (of 
e 5 / 2 £ order) were dropped, as they would be of the same 
order in e as the PN correction of the spin-orbit terms 
[l8j . not included in our treatment. 
Finally we present our concluding remarks in Sec. [V] 



II. LIMITS OF VALIDITY 

During the inspiral, the PN parameter increases. 
Therefore the condition £ < £ x = 0.1 cannot be obeyed 
at distances larger than a certain ri, which depends on 
the particular mass ratio. The value of the correspond- 
ing limiting PN parameter is ei = Gm/c 2 r\ = lOOv 2 . 
Starting from these values of the separation and PN pa- 
rameter, the approximation holds until the PN expansion 
breaks down. Levin, McWilliams, and Contreras argue 
in Ref. (39[, that this should be at £2 = 0.1, and we 
adopt this upper limit. The argument goes as follows. 
At about £2 = 0.1 (well above the innermost marginally 
stable circular orbit) the PN expansion breaks down, as 
the 3.5 PN dissipation term (which is positive) becomes 
comparable to the 2.5 PN order contribution (which is 
negative). Summing up, the regime we study lies be- 
tween £1 and £2. With increasing mass ratio this interval 
shrinks, and it vanishes at f m ax = 0.0316 ks 1 : 32. 

For a specific gravitational wave detecting instrument, 
the lower frequency limit of the best sensitivity band and 
£2 determines the highest possible total mass of a binary 
source for the waves to be detected by the respective 
instrument. Kepler's third law gives a leading order esti- 
mate to the total mass m of the system on circular orbit 
in terms of the gravitational wave frequency / and PN 
parameter e as 

- = ^ 3/2 r 1 . (4) 

The lower frequency bounds are 10 Hz for Advanced 
LIGO [H|, 1 Hz for the Einstein Telescope (HMO" 3 Hz 
for LAGRANGE [H] and 10~ 5 Hz for LISA (H. The 
respective values for m max are thus 202Af© (Advanced 
LIGO), 2O2OM (Einstein Telescope), 2 x 1O 6 M (LA- 
GRANGE) and 2 x 10 s M (LISA). 

Further, a lower limit ^ m i n arises if we fix the smaller 
mass to be that of a neutron star. This leads to ^ m i n ~ 
0.007 w 1 : 143 for Advanced LIGO, v mia 7 x 10~ 4 
for Einstein Telescope, v min « 7 x 10~ 7 for LAGRANGE 
and v min w 7 x lO" 9 for LISA. 

Next we estimate the time-interval spent by the binary 
evolving (due to gravitational radiation dissipation, con- 
sidered to leading order) on quasicircular orbits from £1 
to £2, as function of the mass ratio v e [^min^max]- For 
this we rewrite the leading order radiative orbital angular 
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FIG. 1: (Color online) The color code and contour lines represent the time interval during which a system evolves from 
max (£i , £ f min ) to ei for Advanced LIGO (top left), Einstein Telescope (top right), LAGRANGE (bottom left), and LISA 
(bottom right) gravitational wave detectors, represented as function of the total mass m and mass ratio v. The color code and 
contour lines are displayed logarithmically. The upper bound for the total mass emerges from the lower frequency limit of the 
best sensitivity range for each detector. The (total mass dependent) minimal mass ratio (visible on the top panels) arises by 
assuming the smaller mass to represent a neutron star with typical mass of 1.4 Mq. The maximal mass ratio (the cutoff on the 
right) is determined by the requirement, that E\ = 100^ 2 < 0.1. In the larger part of the allowed parameter space (the region 
on the left of each dashed curve, representing e\ = £/ m in) the gravitational waves entering into the best sensitivity range of the 
respective detectors are described by SDW. In the regions on the right of the dashed curves however SDW can not be applied 
from the lower frequency bound of the best sensitivity range, hence the time intervals are shorter. 



frequency evolution given in Ref. [38[ as 

<t c -t), (5) 



(1 + ^) 2 _ 8/3= 2 8 (Gm) 5/3 



5c 5 

(where t c — t is the time left until the final coalescence) . 
Rewriting Eq. (0| as 



(GmLu) 2/3 



where u) — nf, we obtain 



v bGm 



(6) 



(7) 



Then the time elapsed during the evolution from e\ to £2 
can be calculated as 



At 



5Gm (1 + vf , _ 4 



2 8 c 3 v 



(e l 4 - e 2 



(8) 



Thinking of a specific instrument we can investigate 
the time necessary for the evolution starting from the 



lower frequency range (/min) °f a given detector (when 
this belongs to a value £/ m i n > £1) until £ 2 . On Fig. Q] 
we represent this time interval necessary for the evolu- 
tion from max (£1, e/min) to £2 for the enlisted detectors 
(Advanced LIGO, Einstein Telescope, LAGRANGE, and 
LISA). 



III. SPIN-DOMINATED GRAVITATIONAL 
WAVEFORMS 



Gravitational waveforms to 1.5PN order 



Gravitational waveforms including spin-orbit effects 
were previously calculated in Ref. [6] to 1.5 PN order 
for nearly circular orbits, based on the radiative multi- 
pole moments of 1.5 PN order accuracy, and reproduced 
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(with some misprints corrected) in Ref. [TiJ as 



2Gfi (Gm 



Q'r 



PQ l c J e 



i> L: >o:> + P lb Q% 



Tail 



-3/2 



(9) 



TT 



where TT denotes the transverse trace-free projection 
into the plane orthogonal to the direction N pointing 
from the source to the observer (the line of sight), and 



Q l i = 2 [XX J - ? 



(10) 



pO.SQij = " (6(N.f)f( l A j) 

1 + 1/ L 

+ (N • A) [f i f j - 2X i X j ] | , (11) 



PQl J = 



^(l-3?7)|(N-f) 2 [5Pf 
-16(N-n)(N- X)f (i X ]) 
+ (N ■ A) 2 [3A 4 A J 



7X l X 3 ] 



+-(19-37?)^ - X i X j ) 
3 



2c 1 + v 
Gm 2 l — i/ 



f (i ((a- — S) x N) 



f) 



(12) 



P 15 Ql j = 



1-v 

1 



(1 - 277) 



-(N-A) 3 (f l F-4A l A J ) 



(N-f) 2 (N-A) (58X l X J -37f 



fit; 

- — (N ■ f) 3 f < l A j) + 15(N • r)(N • A) V 4 A J) 



-(N-A) 



^(101 - \2rj)f l f j - ~(19 - 4r7)A l A J 



i(149-36r;)(N-f)f (4 A j) 



2c 



2A (l [fx(o-)] 
-2(N • A) f(ff)xN 



A l A J 



L N -(2S + 3<r) 



L n -(S + <t) 



Ax(4S + 5<r) 



j) 



+2(N-f) [(<t)xn! X j) \ , 



-41n[ — )f<*A j) 



(14) 



Here r % and A 1 are the components of the separation 
vector and of A = Ln x r = v/ru, the parameter i] = 
H/m = v I (1 + v) is the symmetric mass ratio, and 
we introduced the spin combinations S = Si + S2 and 
a = vSi + v~ 1 S2- The last term, P l b Q % ^ aa arising from 
the contributions of the gravitational wave tails was first 
presented in this form in Ref. (H)]. The parameter ljq is 
an arbitrary constant frequency scale [401 ] . Unit vectors 
carry a hat. 

The two gravitational wave polarization states h + and 
h x emerge as the linear combinations 



(15) 
(16) 



of the h iJ tensor components perpendicular to N 6] , with 
y = J x N/ JxN and x = y x N 3 . 

B. The spin-dominated regime 

The total angular momentum J = Ln + Lpn + Lso + 
Si + S2 is conserved to 2PN accuracy [5|. Here Lpn = 
epatLn and Lso are the PN and spin-orbit corrections 
in the orbital angular momentum vector. Starting from 
their expressions Eq. (39) in Ref. [1| and Eq. (B23) 
in Ref. |2l| . and the definition of the PN parameter e we 
find 



^PN 

Ln 
Lso 
Ln 



epN ~ £ 



7 1 £,e 



1/2 



2 2(l + £e 1 / 2 ) 2 

4rv i/2 



Xi 



4(1 + ^/2 

+3 (Xi + X 2 ) + 4^ 1/2 X 2 

Expanding these to linear order in £, we have 

Lpn 
Ln 
Lso 
Ln 



(17) 



(18) 



O (e 2 exi) 



(19) 



(13) 



3 We mention that the polarization vectors x and y are rotated by 
tt/2 in Ref. [T3 |. resulting a global sign difference in the definition 
of h+ and h x . 
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Keeping the terms to e, e£, e 
set of equations, J becomes 



3/2 



orders only in the above 



J= 1 



■Si . 



(20) 



Therefore to the accuracy we are interested in, both the 
second spin and the spin-orbit contribution to the orbital 
angular momentum could be dropped. Then in the tri- 
angle with sides Si, (1 + 7e/2) Ljv and J, containing the 

angles a — arccos • Ln) and /8 1 = arccos • Sx^ the 

sine theorem gives 




D.015 0.02 0.025 0.03 

V 



sin 0! 



7 \ L 



N 



Si 



sin a 



1 j. 7 \ z ■ 

1 H — e — sin a 
2 J XT 



(21) 



As a has no preferred value and £ is small, we conclude 
that Pi is of order £. (From the assumption £ < 0.1 it 
follows that sin/3]^ < 0.1.) Thus any trigonometric func- 
tion of Pi can be approximated by its Taylor expanded 
form, to first order accuracy (dropping £ 2 < 0.01 terms 
and higher). 

It was shown in Ref. [9|, that the angle K\ — a + P l 
stays constant during the inspiral, when the smaller spin 
is negligible. Therefore we replace a everywhere by Ki — 
Pi and then expand to first order in P 1 . We will show in 
the next subsection that the double expansion of h+ and 
h x in e and £ leads to the following structure of terms: 
1, Pi, e 1/2 , e 1/2 Pi, e, e£, e 3/2 . (We note here that the 
leading order contribution from the smaller spin £2 would 
appear only at £ 2 £ 2 orders.) Since e increases during the 
inspiral, the terms Pi, e x l 2 Pi, e£, all kept, increase as 
well, however the e 2 term, discarded, increases at an even 
faster rate. Therefore when either of the enlisted terms 
becomes comparable to the e 2 contribution, we have to 
drop it as well. Figure [2] shows the terms to be kept in 
the parameter space {s,v). 



C. SDW inspiral waveforms 



We give the gravitational wave polarizations {15} and 
(fTTJf in the source coordinate system defined as follows. 
The z,s-axis lies along the direction of J, and the xs~ 
axis is assigned by the projection of N perpendicular to 
J. The polar angle of N will be denoted by 0. In this 
coordinate system the components of the vectors r, A, 
and Si are 



FIG. 2: (Color online) The figure represents the parameter 
ranges v and e for which all terms given in the waveform are 
larger than the pure 2PN (e 2 ) contributions (area 1). In the 
parameter range 2 the contributions e£, eP 1 are becoming 
comparable with the pure 2PN terms, falling outside of the 
limits of validity of our approach, hence they can be neglected. 
In the parameter range 3 all mixed terms in the double expan- 
sion become negligibly small, while in the parameter range 4 
all terms containing £ and /3 1 can be dropped. In the white 
area, the condition £ < 0.1 doesn't hold during the inspiral. 
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= — sin 



3tt 



cos < 



cos a cos 
'3tt 



3tt 



sin i 



— = cos 
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cos a sin 
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3tt 

T 



(22) 
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(23) 



G 



Six = —mj]v Xisin/3 1 cos(0 1 - 0„) , 
c 

G 

Siy = — TO 2 ?7J/ _1 x 1 sin^ 1 sin((/) 1 - 4> n ) , 

Si z = — m 2 T)v~ 1 Xi cos/?! . 
c 



(24) 
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Here the angle (j> 1 is the azimuthal angle of the domi- 
nant spin. (f) n is defined as the angle between the inter- 
section of the planes perpendicular to J and Ln and an 
arbitrary inertial axis taken in the plane perpendicular 
to J, and 4> is the orbital phase. These angles are shown 
on Figs. 1 and 2 of Ref. |20f . The time evolution of these 
angles was given in Ref. [2l| . 

We substitute r, A, and Si from Eqs. §FI§ -(JM]) into the 
expression of h + and h x given by Eqs. (TT5|) and (fTB) . The 
1.5 PN tail contribution contains terms with ln(u;/wo)> 
which can be absorbed into the leading order term by 
redefining the phase (j) as ip — <i> ~ 2e 3 / 2 ln (ui/loq) called 
shifted phase [4(|, pjj. The change of 4> into ip in the 
higher order terms gives modification to the waveforms, 
however these are of higher order than e 3 / 2 . The double 
expansion, including the redshift dependence, of h + and 
h x in e and £ leads to the following structure to £ and 
e 3 / 2 orders: 



2G 2 mV/ 2 £ 
c 4 Dr 



00 



+E 1 ' 2 (h^ + ^h^ -2£h° + 
+e (hi - 4£h° + 5 



h\ so 



+£ 3/2 / fc 1.8 + h 15SO + h lMai> 



(25) 



The indices of the various contributions to the wave- 
form refer to the respective order in the double expansion 
(e,£), respectively f3 1 = 0(£). Thus the terms h+ in- 



elude leading order (e , £ ) contributions, while the h 

X 

terms are the (e ,^) order contributions. The rest of 
the terms without explicit spin magnitude dependence 
are h°+ 5 , of (e°- 5 ,£°) orders, h + 5p , of (e 5 ,^) orders, 

X X 

h\, of (e 1 ,^ ) orders, h^f , of (e 1 ,^) orders and h+ 5 , of 

(e 1 ' 5 ,^ ) orders respectively. The leading order contri- 
bution including the spin mag nitude is h\ so of 

X 

orders, and its corrections h\ and h+ 5SO of (e 1 ,/^) 

X X 

and (e 1 ' 5 ,^ ) orders, respectively. The leading order tail 
term is given by h\ btal \ which is of (e 1 ' 5 ,^ ) orders. All 
these explicit expressions are given in the Appendix. 



by self-spin, quadrupole-monopole 4 and magnetic dipole- 
magnetic dipole contributions in [381 ) . Higher order spin 
contributions were discussed in Refs. [Hj], (T^ |. 

In order to employ these results for the SDWaveform, 
we need to derive the ui (e) dependence to e 2 orders ac- 
curacy. 

Eq. (3) of Ref. [H| (leaving out the magnetic dipole- 
magnetic dipole contributions, which are only important 
for magnetar binaries) gives: 



; 3/2 c 3 

Gm 
^171 



xi 



i 

2 
9 



Xi COS Ki I £ 



-3/2 



- — + — COS Kl 



(26) 



Note that the spin-spin terms are shifted to e 5 / 2 £ order, 
thus are neglected. 

The radiative orbital angular velocity evolution (w) is 
given by Eq. (7) of Ref. [HI]. This evolution for un- 
equal mass case at linear accuracy in £ and e 2 (with the 
neglection of e 2 £) becomes 



96 e 6 £c 6 
5" (Gmf 



+ 4tt 



79 



12 
697 465 



1105 

~336~ £ 

23 



-Xi cos ki I e 



-3/2 



9072 



-xl 



35 



325 . 2 

sin Ki 

96 16 



(27) 



Integrating Eq. (7) of Ref. [Hi and employing Eq. (|2"6"|) 
we find to linear accuracy in £ and to e 2 orders beyond 
leading order: 



e 

256 



265 ^ 
252 £ 



62 



"Xi COS Kl 



^ 32 

5£ - — 7T + - 

5 

24 804 463 
508 032 



63 577 . 2 

+Xi sm Ki 

48 



.3/2 



e 



Here r is a dimensionless time variable defined as 

„3 



T : 



5G 



m 



1 - 2e 



1/2 e) (*c - t) 



(28) 



(29) 



IV. THE PHASE OF THE GRAVITATIONAL 
WAVEFORM IN THE SPIN DOMINATED 
REGIME 



where t c is the phase and the time at the final coales- 
cence. Integrating twice Eq. (7) of Ref. [38[ gives the 



The radiative evolution of the orbital angular veloc- 
ity up to 2PN orders.including spin-orbit and spin-spin 
effects was given in [6j. Later, this was complemented 



4 The coefficient p_i of Ref. 38] is pi = — G 2 /e 4 aJx 2 m 2 /m 2 , where 
uj = 1 for black holes. 
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accumulated orbital phase as 

^ 3 < 1195 



= L-/l + 2 E V2£. 



1008 



-e 



3925,. 175 ^ 3/ , 

~ + "504~ C + ~T Xl C ° S Kl ' 

21440 675 
~ 1016 064 



., ,'375 3425 . , 
+Xl (^-^6- Sm Kl 



where C is the phase at the final coalescence. 



(30) 



V. CONCLUDING REMARKS 

In this paper we have analyzed spinning compact bi- 
nary inspiral in the case when one of the binary compo- 
nents has much larger mass than the other one and it 
spins fast. In such cases the end of the inspiral is char- 
acterized by a much larger spin of the dominant binary 
component as compared to the orbital angular momen- 
tum, such that the parameter £ = Xi {Ln / S\) < 0.1. 
This latter condition selects the mass ratios 1 : 30 and 
below. Similar considerations lead to the conclusion, that 
the second spin can be neglected. Therefore, from all 
available angular momenta of the system, the larger spin 
dominates. Hence we call the corresponding waveforms 
spin-dominated waveforms (SDW). 

The SDW were given as a double expansion in the PN 
parameter e and the smallness parameter £. The angle /3 X 
between the dominant spin and total angular momentum 
is also small, of order £ (this is because the larger spin 
gives the main contribution to the total angular momen- 
tum) The usual mass variables employed for gravitational 
waveforms, the chirp mass M c hi rp = (mim2) 3 ^ 5 to -1 / 5 = 
mr] 3 / 5 , and the mass ratio v relate to the total mass m 
and £ as follows: 



1 : 140 for Advanced LIGO, 7 x 10~ 4 for the Einstein 
Telescope, 7 x 10" 7 for the LAGRANGE and 7 x 10~ 9 
for the LISA missions, respectively. 

The SDW amplitudes were presented to (e 3 / 2 ,£) or- 
ders; while the phase of the gravitational waves to (e 2 , £) 
orders (omitting the highest order mixed terms). To this 
accuracy the amplitude includes the (leading order) spin- 
orbit contributions. The phase includes the (leading or- 
der) spin-orbit, the self-spin and the mass quadrupole- 
mass monopole contributions. Note that expressing the 
mass ratio in terms of £ shifts the usual e order, and this 
is what happens with the spin-spin contributions, which 
would appear only at higher orders. This is consistent 
with the smaller body's spin being insignificant in the 
chosen parameter range in comparison with the larger 
spin and orbital angular momentum. 

A comparison showed that the derived SDW, when 
written with the same level of detail as the waveforms 
given in Appendix A of Ref. [l4j , are approximately 80% 
shorter. This is due to the smaller number of variables 
and most importantly, the second series expansion we 
employ. We expect that the SDW will turn useful both in 
modeling gravitational waves emitted by binary systems 
consisting of an astrophysical black hole with a smaller 
black hole or a neutron star companion; and for super- 
massive black hole encounters. 
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chirp 



3/5 



me 



3/10^3/5 
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The limits of validity of these waveforms were studied, 
their typical lengths in the instrument best sensitivity 
range were estimated for the Advanced LIGO, the Ein- 
stein Telescope, the LAGRANGE and the LISA gravita- 
tional wave detectors, respectively. While the SDW hold 
for any mass ratio smaller than 1 : 30, lower bounds for 
the mass ratios were derived from the best sensitivity 
frequency range expected for these forthcoming gravita- 
tional wave detectors. The mass ratio lower bounds are 



Appendix A: Explicit Spin-dominated waveforms 

We give below the various contributions to the wave- 
form expressed in the double expansion (e, £), remem- 
bering that Pi = O (£). The (e°, £° ) contributions are 



4/i° 



2h u 



J2 [(sin 2 9-2) cos(20„ ± 2^)c[ ±0) 

—2 sinKi sin 26>sin(0„ ±2V>)fc (±) 

+6 sin 2 ki sin 2 9 cos 2ip , (Al) 

J2 [cos 9 sm(2(t) n ±2iJj)c[ ±0) 

-2sin0sin/cicos(0 n ± 2^)fc (±) j , (A2) 
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while the (e°, (3{) contributions read 



2hf SO 



h 



opso _ 



J2 [sin20sin(0„±2?/;)4 ±O) 
+,- L 

+ sin ni (sin 2 0-2) cos(20„ ± 20)fc (±) 
— 3 sin 2ki sin 2 9 cos 2-0 , (A3) 
J2 [cos0sin/tisin(20 n ±2V>)fc (±) 



TABLE I: The terms <4 ±n) in all PN orders, n is the PN 
order, i is the number of the constant in the appropriate PN 
order. 



+ sin 9 cos(0„ ± 2V>)4 ± °' ) ■ 



(A4) 



The (e 5 , £°) corrections to the gravitational wave- 
forms h+ are 



4 cos 9 sin ki sin 2 [—45 sin 2 ki cos 3V> 

+ cos -0(5 sin 2 Ki — 4) ] 

+ { — sm2 K i sin [(sin 2 9 



-2) (sin(30„ ± 0)fc< : 
-45 (2-3 sin 2 6) (sin(0„ ± 30)fc (±) 



+ sin0 



(sin(0„±0)ci ±O - 5) ) 



-9 (sin 2 6>-2) (sin(30„ ± 30)4 ±O ' 5) ) 
+2 cos (9 sin ki (cos(20„ ± V>)4 ±0 ' 5) ) + 9 ( 2 
-3sin 2 6») (cos(20„ ±30)4 ±O) ) }, (A5) 



32/i 



0.5 



— 16 sin 2ki sin -0 sin 2 

+ {cos sin 2 m sin9 [cos(30„ ± ip)k^ 
+,- 

+45 (cos(0„ ±3^)fc (±) ) 

+ isin20[(cos(0 n ±^)4 ±o - 5) ) 

+9 (cos(30„ ±3^)4 ±0 ' 5) ) 



+2sinKi 



(sin(20„±0)4 ±O - 5) ) 



-9cos2(9(sin(20„±3VO4 ±O) )]} » ( A6 ) 



T 2fe (±) 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 



sin 2 0) 



0.5 1 4fc (±) (6 

4fc<±> 
±2 (6 - sin 2 0) fc (±) 

12fc(±) 
±2fc( ± )(2 sin 2 -3) 

-2fc<±) 
T 24 ±0) (6 - sin 2 0) 
±2Jfe< ± > 

34 sin 2 ± 2 (5 sin 2 - 46) cos ki 
22 ± 46 cos ki 

-2fc (±) (3 -2 sin 2 g) 

±8fc<±) 
6fc (±) (sin 2 + 5) 



44 - 



sin 2 0) 



2fe(±) (4 - 
±2fc (±) (2 sin 4 6» + ll sin 2 - 38) 
6fc (±) (3 sin 2 + 5) 
±2fc (±) (4sin 2 + 19) 
-2fc (±) (3 sin 2 6» -4) 
±2fc (±) (4 - sin 2 0) 
±6fc (±) (5 + sin 2 9) 
T 4fc<±) 

fc (±) (22 + 29 sin 2 - 16 sin 4 0) 
2fc<±) 
±6fc (±) (3 sin 2 6» + 5) 
-fc (±) (20 sin 2 + 11) 
=p2fc (±) (3 sin 2 9 - 4) 



1.5 1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 



±4fc (±) (sin 2 9-6) 
±4fc< ± > (sin 4 9 + 42 sin 2 0-166) 
16/fcW 

8fc (±) (sin 4 + 8 sin 2 - 28) 
8fe (±) (-332 + 94 sin 2 + sin 4 0) 
±8fc (±) (38 - 42 sin 2 0-9 sin 4 0) 
-16fc (±) (l52 - 46sin 2 0-9sin 4 0) 

±24fc (±) (3 sin 2 - 10) 
-8fc (±) (160 - 204 sin 2 - 63 sin 4 0) 
±4fc (±) (3 -2 sin 2 0) 
-8fc (±) (14 + 3 sin 2 0) 
-16fc (±) (15 sin 2 + 76) 
-8fc (±) (5 sin 2 + 166) 
-8fc (±) (80 + 63 sin 2 0) 
±4fc (±) (166 - 125 sin 2 0-8 sin 4 0) 
+8fc (±) (38 - 61 sin 2 - 24 sin 4 0) 
±8fc (±) (5 -4 sin 2 0) 



while the (e , /3 : ) contributions read 
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TABLE II: The terms bg" J in all PN orders, n is the PN 
order, i is the number of the constant in the appropriate PN 
order. 



n 


i 


A±n) 
°(i)0 





1 


-1 




2 


-2 


0.5 


1 


± (46 - 5 sin 2 0) 




2 


±3 




3 


2-3sin 2 




4 


±23 




5 


- cos 26» 




6 


±4 




7 
8 


— 
+3 




9 


-15 (2 -3 sin 2 6 1 ) 




10 


-15 




11 


+4 (3 - 2 sin 2 0) 


1 


1 


8 




2 


+ 18 ±7 sin 2 




3 


+3sin 2 0±6 




4 


-22 - 29 sin 2 + 16 sin 4 




5 


±26 sin 2 + 18 




6 


±11 ±20 sin 2 




7 

8 


+6sin 2 0±6 
2 (11 -5 sin 2 6») 




9 


6 (7 ±9 sin 2 6») 




10 


-11 




11 


+3 (8 -20 sin 2 9 ±7 sin 4 9) 




12 


±3 




13 


3 (19 sin 2 6 ± 14) 




14 


±12 cos 26» 




15 


22 -21 sin 2 9 


1.5 


1 


5sin 2 0-6 




2 


18 sin 4 6 ±252 sin 2 9- 188 




3 


±20 




4 


±9 sin 4 9 + 90 sin 2 9 ± 56 




5 


+4 (1184 - 172 sin 2 - 7 sin 4 0) 




6 


2 (46 ±48 sin 2 0-99 sin 4 0) 




7 
8 


+ 12 (72 ± 110 sin 2 (9-63 sin 4 6) 
144 (sin 2 0-2) 




9 


+3 (-204 ± 406 sin 2 9 - 189 sin 4 9) 




10 


3-4 sin 2 9 




11 


±28 + 31 sin 2 9 




12 


+432 + 876 sin 2 9 




13 


+4 (71 sin 2 6 ± 592) 




14 


±306 + 651 sin 2 9 




15 


2 (94 - 173 sin 2 0- 24 sin 4 9) 




16 


-2 (46 ± 25 sin 2 9 - 180 sin 4 0) 




17 


48 cos 2 6 



TABLE III: The terms dffi in 0.5, 1 and 1.5 PN orders, 
n is the PN order, i is the number of the constant in the 
appropriate PN order. 



n 


i 


d (±n) 


0.5 


1 


5 (3 sin 2 9 -2) 




2 


-1 




3 






4 


-5 




5 






6 


3 




7 
8 


-3 (2 -3 sin 2 9) 




9 
10 


_ 




11 


3 cos 26» 


1 


1 


±4 




2 


6-14 sin 2 9 




3 


2 (sin 2 6»-l) 




4 


q=2 (8 - 20 sin 2 9 + 7 sin 4 6») 




5 


6 - 7sin 2 9 




6 


=Fl6sin 2 6»±8 




7 


3 sin 2 9 - 2 




8 


±9 (sin 2 0-2) 




9 


±3 (18 -7 sin 2 9) 




10 


±9 




11 


4(2-8sin 2 6»±7sin 4 9) 




12 


-2 




13 


±6 (9- 13 sin 2 6») 




14 


2 (7 sin 2 9 - 2) 




15 


Tl8cos 2 6> 


1.5 


1 






2 


+6 sin 4 9 ± 72 sin 2 + 20 




3 


-12 




4 


-15 sin 4 ± 22 sin 2 0-8 




5 


1920 - 2832 sin 2 - 84 sin 4 




6 


±6 (10 - 44 sin 2 ± 27 sin 4 0) 




7 


-4 (88 - 422 sin 2 ± 171 sin 4 0) 




8 


±12 (14 -9 sin 2 0) 




9 


-9 (28 - 126 sin 2 ± 105 sin 4 0) 




10 






11 


9 sin 2 0-4 




12 


-176 ± 756 sin 2 




13 


12 (7sin 2 0±8O) 




14 


-126 ± 189 sin 2 




15 


±2 (10 - 41 sin 2 ± 36 sin 4 0) 




16 


+6 (10 - 49 sin 2 ± 44 sin 4 0) 




17 


+4 (7 -8 sin 2 0) 
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TABLE IV: The terms in 1 and 1.5 PN orders, n is the 

PN order, i is the number of the constant in the appropriate 
PN order. 



6 3 (6 - 36 sin 2 + 45 sin 4 0) 

7 +3 (232 - 510 sin 2 + 243 sin 4 9) 

8 9 (6 -5 sin 2 6) 
9 

10 
11 
12 
13 
14 
15 
16 
17 



+348 ± 591 sin 2 9 
± (273 sin 2 9 - 118) 

(2 - 13 sin 2 6» + 12 sin 4 0) 
-9 (2 - 13 sin 2 9 + 12 sin 4 9) 
-3 (3 -4 sin 2 9) 



TABLE V: The terms dffi in 1.5 PN order, n is the PN 
order, i is the number of the constant in the appropriate PN 
order. 



n 


i 


fe (±n) 

Ml 


n i 






1 


1 


-1 


1.5 1 








2 


_ 


2 




— 




3 


_ 


3 




1 




4 


-2(2-8sin 2 + 7sin 4 9) 


4 




— 




5 


_ 


5 


28- 


126 sin 2 + 105 sin 4 




6 


2 - 7sin 2 (9 


6 




— 




7 
8 


— 

-8 cos 2 9 


7 

8 


27 (8 


-22 sin 2 + 15 sin 4 0) 

— 




9 


8 (3 - 7 sin 2 9) 


9 




— 




10 


4 


10 




- 




11 


_ 


11 




— 




12 


- 


12 




-243 sin 2 + 108 




13 


-4 (7sin 2 <9-6) 


13 




7 (2 -3 sin 2 0) 




14 




14 








15 


-4 (2 -3 sin 2 9) 


15 






1.5 


1 




16 








2 


-15sin 4 0+12sin 2 0-2 


17 








3 


+5 










4 












5 


+ (236 - 294 sin 2 + 21 sin 4 0) 









64h° + 5 ' = 4 cos ki cos 9 sin 2 9 [135 cos (3V>) sin 2 «i 

+ cosV>(4- 15sin 2 Ki)] + £ {2 cos 9 [9 (2 

+,- 

-3 sin 2 9) (cos(20„ ± 3^)4 ±0 - 5) ) 



+ cos(2</>„ ± ^)c,." 



(±0.5) 



sin m sin# (sin 2 9 



-2) 



27 



sin(3^±3V-)4 ±0) ) 



• sin(3<£„ ± V)4 : 



±0.5) 



sin Ki sin9 [45 (2 



-3 sin 2 9) (sm{4> n ± 3V>)4 ±0 ' 5) ) 

+ sin(0„±^4 ±O - 5) ]} , (A7) 



32/i°' 5/3 = 32sin-0sin 2 6»cos2Ki 



+ 



cos sin K\ sin 



27cos(30„ ±3-0)4 ±o) 



+ [cos(3</>„ ± V) + 45 cos(0„ ± 3^)] 4 ±0 ' 5) 

+4o °' 5) cos(0„ ± V) 

- [18 cos 29 (sin(20„ ± 3ip)c%) 

-2sin(20„±^)4f- 5) ] } . (A8) 

The (e 1 , £ ) corrections to the gravitational waveforms 
/i+ are 
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48h\ = 2sin 2 6»sin 2 Ki [sin 2 Ki5 (7sin 2 6» 
—6) (cos 2-0 — 4 cos 40) 
-2 (15 sin 2 9 + 51) cos 20] 
+ { 16 sin 3 ki sin 26* (7 sin 2 



+,- 



-3) sin(0„ ± 40)fc (±) - (sin 2 9 
-2) sin 2 sin 2 9 cos(40„ ± 20)4 ±O) 
+4 (sin 2 9-2) sin 2 9 cos(40„ ± 40)4 ±1} 
+2 sin Ki sin 29 sin(0 n ± 20)4 ±1 ^ 

+ sin(30„±20)4 ±1) 

-8 cos 2 9 sin(30„ ± 40)4 ±O ' 5) 

+2cos(20„±20)4 ±1) 

- 16 sin 2 m [7 sin 4 9-2(4 sin 2 9 

-l)]cos(20 rl ±40)4 ±o) } , 



(A9) 



24hi 



60 cos Ki cos 9 sin 2 ki sin 20 sin 2 9 

+ [8 sin 3 /ti sin 6> (7 sin 2 9 

+,- 

-6) cos(<?!>„ ± 40)/c (±) 

— cos sin 2 ki sin 2 sin(4</>„ ± 20)4 ± °' ) 

+2 sin Ki sin 9 cos(</>„ ± 20)4 ±1 ' ) 

+4 cos 6» sin 2 sin(40„ ± 40)c^ ±:L) 

+2cos6»sin(20„ ±20)4 ±1} 

-8 cos 9 sin 2 ki (7 sin 2 9-2) sin(20„ ± 40)4 ±O) 

-8 (2 - 3 sin 2 (?) sin /ti sin 9 cos(3</>„ ± 40)4 ±O - 5) 

+2sinKisin0cos(3(?!)„ ±20)4 ±1} j , (A10) 



2Ah^ = cos Ki sin ki sin 2 #[4 (15 sin 2 6> + 51) cos 20 
+20 sin 2 ki (7 sin 2 - 6) (4 cos 40 - cos 20)] 

+ J2 {sin20sin(30„ ±20)4 ±1} 

+ sin20sin(0„ ±20)4 ±1} 

+8 cos 2 9 sin 29 sin(3</>„ ± 40)c^ 1} 

+2 sin ki cos(2^„ ± 20)c^ 1} 
+8sin 2 Kisin26>(7sin 2 6> 

-3) sin(0„ ± 40) [±3fc (±) + 4 sin 2 k x 

+16shiKi (2 - 8sin 2 6» 

+7 sin 4 9) cos(2</>„ ± 40)c^ 1} 

+ sin ki sin 2 9 (sin 2 9-2) cos(40„ ± 20)c^ 1} 

~8cos(4,/)„±40)4 ±o - 5) ]} , (All) 

12h]f = 30 cos 9 sin ki sin 20 sin 2 9 (3 sin 2 k\ 

-2)+ £ {sin0cos(0 n ±20)4f 1) 
+,- L 

+2 cos sin K\ j^sin(20„ ± 20)c^ 1 ^ 
+4 (7 sin 2 6>-2) sin(20„ ± 40)4f 1} 
+ cos 9 sin ki sin 2 9 sin(4</>„ ± 20)42 ^ 
-8sin(40„±40)4 ±o - 5) " 
+ sin<9 cos(30„ ± 20)4g 1} 
+4 (2 - 3 sin 2 9) cos(30„ ± 40)4 1} ] 
+4 sin 2 Kl sin (6 - 7 sin 2 9) (r3fc (±) 
-4sin 2 Kl ) cos^ ±40)} . (A12) 

The (e 1 , £° ) and (e 1 , £) SO corrections to the gravita- 
tional waveforms h+ are 

X 

h\ so = ^ sin 9 Y. [±fc (±) sin(^ ± 0)] , (A13) 
1 +,- 



ISO 



Xi 



sin {2 sin Ki sin sin 



+ cos(9£ [±fc(±) cos(</>„ ± 0)] > ,(A14) 



while the (e 1 , (3^ contributions read 



— h 1 / 30 = V {t^cos 9 sm(2<p n ±ip-(j> 1 ) 
Xi +,- 

+ sin Ki sin 6* [± sin(0„ ± 0) 
±cos(0„±V-^i)]} • (A15) 
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Xi 



ft 1 / 50 = - [2 (sin ^ + 2) cos «i sin ^ 

+2 cos cos 0] sin 2 

+ ]C {2cos#sinKi sin0 [±cos(0 n ± 0) 



,(±i-5) 



+,- 



+ sin(0„±0-</> 1 )]±(sin 2 
-2)/c ± cos(2^±^-0 1 )} , 



(A16) 



6144/i^ 5 = 192cosKisinKisin 2 6'[sin^(64 
- sin 2 ki (7 sin 2 - 6) + 4 sin 2 9] 
+27 sin 3^ sin 2 ki (7 sin 2 0-6)] 

+ J2 ( + 4 si n 3 «i sin 2 sin(40„ ± 0)4^ 
+,- L 

—2 cos sin 4 ki sin 3 cos(50„ ± tp)k^ 
-243 cos sin 2 m sin 3 cos(50„ ± 3-0)4 ±O - 5) 
-625 cos sin 3 cos(50„ ± 5^)4 ±1-5) 
+3 sin 2 /ti sin20cos(30„ ± f/O c n * 5 ' 
+27cos0sin0cos(3</>„ ± 3 , 0)c^ 1 ' 5 ^ 



The (e , £ ) corrections to the gravitational waveforms 
h+ are 



+1875 cos sin 2 «i sin (4 
-9 sin 2 0) cos(30„ ± 50)4^ 
+2sin20cos(0„ ± V)4f 1 ' 5) 



12288/^ = 12 cos sin ki sin {cos Si/ 1 [1701 (2 

-3 sin 2 6») sin 4 «i + 72 sin 2 m (63 sin 2 9 
+178)] + cos0 [-14 (2 - 3 sin 2 0) sin 4 K\ 
—8 sin 2 (7 sin 2 9 + 162) + 16 (sin 2 + 66)] 

-4375 (2 - 3 sin 2 6) sin 4 «i cos 50} + ^ [2 (sin 2 

+,- 

-2) sin 4 ki sin 3 0fc (±) sin(50„ ± V) 
+4 cos 9 sin 3 m sin 2 9 cos(40„ ± V)4 ±1-5) 
+16 cos 6* sin «i cos(20„ ± ip)^ 1 ' 5 ^ 
+1250 sin 4 /ci sin 9 (105 sin 4 

-126 sin 2 + 28) fc (±) sin(0„ ± 50) 

+625 (sin 2 0-2) sin 3 9 sin(50„ ± 5^)4 ±1 ' 5) 

+6 sin 2 «i sin sin(30 n ± ^)4 ±1-5 ^ 

+243 (sin 2 0-2) sin 2 m sin 3 9 sin(50„ ± 30)4 ±O ' 5) 

+4sin0sin(0 n ±f/')4 ±1 ' 5 ^ 

+5000 cos sin 3 m (15 sin 4 - 12 sin 2 

+2) cos(20„ ± 50)4 ±O) 

— 1250 cos sin ki sin 2 9 (5 sin 2 



+27 sin 2 K\ sin 20 cos(0„ ± 3ip) c iA 1 ^ 
+8 sin Kl sin(20„ ± 0)4f 1 ' 5) + 4375 (2 
-3 sin 2 0) sin 4 m sin 20/cW cos(0„ ± 5^) 

+108 sin ki sin(20„ ± 3^)4e 1-5) 
+162 sin m sin 2 sin(40„ ± 30)4y 
-2500 sin 3 ki (2 - 13 sin 2 

+12 sin 4 0) sin(20„ ± 5^)c[ ±0) 
- 1250 sin ki sin 2 0(3 

-4sin 2 0)sin(40„±50)4 ±1) } i ( A18 ) 

The (e 1 ' 5 , £° ) SO corrections to the gravitational wave- 
forms h+ are 



Xi 



-6) cos(40„ ± 5^)cJ ±1} 
+1875 sin 2 K\ sin (8 - 22 sin 

+15 sin 4 0) sin(30„ ± 5^)4 ±0 ' 5) 

+216 cos sin ki cos(20„ ± 3^)4 ±1-5 ^ 
„(±!-5) 



—ft 

Xi 



1.550 



+27sm0sin(30 n ±3^)cf 
+54cos0sinKi sin 2 0cos(40 n ± 3^>)4 ±1-5 ^ 
+54 sin 2 Ki sin sin(0 n ± 3^)4 ±1-5 ^ 



(A17) 



4 sin «4 [cos Ki sin K\ cos 20 n 

— cos 2ki cos sin <f> n sin 
+ cos Ki sin /ti sin 2 (6 sin 2 ip 
-2 + sin 2 n )] 

+ £ {2cos0sinKi sin [(+3*^) 

+,- 

-4 sin 2 «i) sin(0„ ± 2ip)] 
+ cos(20„±20) -2fc (±) + (2coski 
+3) sin 2 ki] (sin 2 - 2) } , (A19) 

—2 cos 4> n sin K\ (sin cos 2/ti 

+ cos sin 2ki sin n ) 

+ {cos0sin(20„ ± 20) [-2fc( ± ) 
+,- 

+ (2 cos K\ + 3) sin 2 «i] 
+ sinKisin0 +3fc' ± - ) 
-4 sin 2 ki] cos(0„ ± 20) } , (A20) 
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The (e 1 ' 5 , £° ) tail corrections to the gravitational wave- 
forms h+ are 



also appearing in the various PN orders of h+ are struc- 

X 

tured as 



, l.dtail 



6 sin Ki sin 8 cos 2vb 



£ cos(2</)„ ± 2^)cp 0) (sin 



-2fc (±) sinKi sin26>sin(0„ ± 2ip)] ,( A21 ) 



±n (±n) . • 2 

W W 



6 (±n) 



j=0 



+d\f ) "> cos j sin 2j Ki . (A23) 



7T 



1.5taz/ 



£ cos6»sin(20„ ± 2vj)c\ 



(±o) 



-2 sin ki sin 6>fc (±) cos(0„ ± 2-0) (A22) 



In the above expressions we have introduced the nota- 
tions fcW = coski =F 1- The constant coefficients c- ± ™' ) , 



Here n denotes the PN order of the coefficients and i, 
7 are serial numbers. The contributions ay^ n , b/^„ , 

djj^, feg^ and dgj° are given in the Tables HI HTl ITTT1 [TVl 
and |V] Note that in the nonprecessing case (ki = 0,tt) 
the contributions (with k + — 0, hr = 2) fully 

define the coefficients c[ ±n \ 
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